Particle decays during the cosmic dark ages 
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We consider particle decays during the cosmic dark ages with two aims: (1) to explain the high 
optical depth reported by WMAP, and (2) to provide new constraints to the parameter space for 
decaying particles. We delineate the decay channels in which most of the decay energy ionizes and 
heats the IGM gas (and thus affects the CMB), and those in which most of the energy is carried 
away — e.g. photons with energies 100 keV < E < 1 TeV — and thus appears as a contribution to dif- 
fuse x-ray and gamma-ray backgrounds. The new constraints to the decay-particle parameters from 
the CMB power spectrum thus complement those from the cosmic X-ray and 7-ray backgrounds. 
Although decaying particles can indeed produce an optical depth consistent with that reported by 
WMAP, in so doing they produce new fluctuations in the CMB temperature/polarization power 
spectra. For decay lifetimes less than the age of the Universe, the induced power spectra generally 
violate current constraints, while the power spectra are usually consistent if the lifetime is longer 
than the age of the Universe. 

PACS numbers: 98.80.-k, 98.70.Vc, 95.30.Cq, 95.35.+d 
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I. INTRODUCTION 

A large correlation between the temperature and E- 
type polarization at large angular scale (low I) was re- 
cently observed bythe Wilkinson Microwave Anisotropy 
Probe (WMAP) [l|. This is a unique signature of re- 
scattering of cosmic microwave background (CMB) pho- 
tons at redshifts relatively low compared with that of 
the last-scattering surface at z « 1100 2]. The required 
optical depth of r e ~ 0.17 can be achieved if reioniza- 
tion occurs at a redshift of z re ~ 20. Although there 
are theoretical uncertainties, such a reionization redshift 
is difficult to reconcile with the star-formation history 
expected in the ACDM model Q, which generally fa- 
vors a reionization redshift of 7-12 Q. Furthermore, the 
thermal history of the intergalactic medium (IGM) con- 
tains further evidence for late completion of reionization 
H . This potential conflict between the evidence for early 
and late reionization might be partially resolved in the 
double-reionization model, where an early generation of 
massive, metal-free stars were formed and partly ionized 
the Universe 0. Nevertheless, even in this model, it is 
not easy to achieve such a high optical depth Q]. 

In light of this, it is worthwhile to consider possible 
alternatives. For example, it has been suggested that a 
high optical depth might be achieved if primordial den- 
sity fluctuations are non-Gaussian 0. Here we consider 
another option. While stellar photons must have con- 
tributed to reionization, it remains possible that other 
energy sources also contribute. Decay of an unstable 
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particle, for example, provides such an alternative en- 
ergy source. In this scenario, a decaying particle, pos- 
sibly part of the dark matter, releases energy during its 
decay, which contributes to the ionization of the IGM. 
Another widely discussed possibility is the radiative de- 
cay of an active neutrino, which might play a role in 
a number of astrophysical phenomena 0, IToj ■ Al- 
though the parameters of the original model are now 
excluded by observations there are still other re- 
gions of decaying-neutrino parameter space, and there is 
no lack of other particle-physics candidates; e.g., unsta- 
ble supersymmetric particles |l2j , cryptons Il3|. R-paritv 
violating gravitinos |l4j , moduli dark matter 1 151. super- 
heavy dark-matter particles 0, 0], axinos fig . ster- 
ile neutrinos |l9j |. weakly interacting massive particles 
(WIMPs) decaying to superweakly interacting massive 
particles(SWIMPs) l20ll. a nd auintessinos |2lj . Recently, 
Hansen and Haiman |l9| suggested sterile-neutrino de- 
cay as a source of reionization. In addition to decaying 
particles, evaporation of primordial black holes [2^ and 
decay of topological defects such as cosmic strings and 
monopoles are also possible source of extra energy input. 
The decay of an unstable particle may also help explain 
the presence of dwarf spheroidal galaxies in the Local 
Group, resolve the cuspy halo problem in ACDM models 
|23l 124 |25| , and serve as a possible source of the ultra 
high energy cosmic rays [2(|. 

From a cosmological perspective, it is particularly in- 
teresting to consider the rich variety of ionization histo- 
ries offered by the particle-decay scenario. In these sce- 
narios, the Universe is not necessarily fully ionized; in- 
stead, particle decay may ionize only a small fraction of 
the gas. If the process lasts for a large range of redshifts, 
it may still contribute a large fraction of the measured 
free-electron optical depth. The presence of a not signifi- 
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cantly damped first acoustic peak in the CMB anisotropy 
spectrum suggests that particle decay should not sig- 
nificantly delay the recombination process at z ~ 1100 
[27L l28l ] , but is it possible that the Universe become par- 
tially ionized during the cosmic "dark ages" at redshifts 
of ten to a few hundred? What is the observational sig- 
nature of such an ionization history? Can this scenario 
be distinguished from late reionization by CMB obser- 
vations? Particle decay may also produce energetic pho- 
tons; can observation of cosmic 7-ray backgrounds place 
constraints on this scenario? 

In this paper we consider these questions. Since at 
low redshift stars and quasars emit ionizing photons, and 
since at the epoch of recombination there is no significant 
increase of entropy, we shall focus mostly on particle de- 
cays in the redshift range between 1000 and 20. Such 
particles produce an optical depth r ~ 0.17 by partially 
rcionizing the Universe at redshifts much higher than the 
value, z ~ 20, required if the Universe becomes fully 
ionized by early star formation. We calculate the CMB 
temperature/polarization power spectra induced by this 
alternative ionization history and show that it can be 
distinguished from the full-rcionization scenario with the 
same r. In some regions of the decay-particle parame- 
ter space, the induced power spectra conflict with those 
observed already, but there are other regions where de- 
caying particles can provide the required optical depth 
and maintain consistency with the measured power spec- 
tra. 

While investigating decaying particles as contributors 
to cosmic reionization, it becomes clear that new CMB 
constraints to the ionization history provide new con- 
straints to the parameter space for decaying particles. 
To a first approximation, the energy injected by particles 
that decay with lifetimes between the ages of the Universe 
at recombination and today either gets absorbed by the 
IGM, or it appears in diffuse radiation backgrounds pflj . 
In the latter case, observed radiation backgrounds have 
traditionally been used to constrain the parameter space 
that consists of the decay-particle lifetime and density as 
well as the energy of the decay products. As we detail 
below, new CMB constraints to the ionization history 
can now provide complementary new constraints to the 
regions of parameter space where the decay energy goes 
to heating and ionizing the IGM. 

This paper is organized as follows: In the next Sec- 
tion, we discuss how energy is dissipated for various de- 
cay channels and what fraction of energy is eventually 
used for ionization, heating the gas, or carried away by 
escaping photons and neutrino. We also discuss how this 
energy is deposited as a function of redshift. In Section 
III, we describe how to calculate the ionization history 
and CMB anisotropy with extra energy input from de- 
caying particles, and we discuss how the result depends 
on the property of the particle. In Section IV, we obtain 
constraints to the decay-particle parameter space from 
the CMB and diffuse backgrounds. We summarize our 
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FIG. 1: Energy loss rate by photons per Hubble time. The 
solid (blue), dotted (green), dashed (red) line are for redshift 
z + 1 = 10, 100, and 316, respectively. We took x e = 0.01 
here, but the results are insensitive to x e for E > 1 keV. 



results in Section V. 



II. DECAYING CHANNELS AND ENERGY 
DISSIPATION 

Depending on the nature of the decaying particle, the 
decay products may include gauge bosons, charged lep- 
tons, neutrinos, quarks, or other more exotic particles. 
These particles may then subsequently decay further into 
other particles, or they may interact with particles in the 
IGM. With sufficient energy, a shower of particles is cre- 
ated. In the end, stable, weakly-interacting particles like 
neutrinos escape, while other particles lose a significant 
part of their energy during the interaction with the pri- 
mordial gas or cosmic microwave background. Some of 
this energy can go into ionizing the IGM, and the effi- 
ciency of converting the decay energy to ionization energy 
is process dependent. Here wc review some of the more 
generic features; in particular, we consider the efficiency 
of converting the rest mass of the decaying particle to 
ionization energy, Xi = Ei/Mxc 2 , where M x is the mass 
of the decaying particle. 



While this paper was being prepared, two preprints 1 3 Ot l.'jll on 
similar questions appeared. Our results agree with theirs where 
our calculations overlap. 
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A. Photons 

In this paper we are mostly interested in the "dark 
ages", 10 < z < 1000, where most of the gas is neutral. 
Photons with energy smaller than 13.6 eV cannot ionize 
hydrogen atoms in the ground state, but if there is a large 
presence of hydrogen atoms in excited states, e.g. at the 
end of the recombination era, z ~ 1000, photons with 
energy E < 13.6 eV may contribute to the ionization. 
When most of the atoms fall to the ground state, photons 
with energy less than 13.6 eV will escape. 

Ultraviolet and soft X-ray photons with energy 
13.6 eV — 1 keV have large photoionization cross sec- 
tions and are largely absorbed locally. For photons with 
energy E > 40 eV, neutral helium absorption dominates 
the absorption, and photoelectrons are produced in the 
process. The photoelectron carries the remaining energy 
of the photon. From here on the energy dissipation pro- 
cess for the initial photon is the same as that for an initial 
energetic electron. 

The absorption processes of hard X-ray and 7-ray pho- 
tons at cosmological distances were discussed in Ref. [32| . 
The processes in which photons can be absorbed or lose 
energy include (i) photoionization of atoms, (ii) Compton 
scattering on electrons, (iii) production of pairs on atoms, 

(iv) production of pairs on free electrons and free nuclei, 

(v) scattering with background photons, and (vi) pair 
production on background photons. In some of these, 
such as Compton scattering, a photon loses only a small 
fraction of its energy in a single scattering event, while 
in others it can lose a significant part of its energy. 

In Fig. ^ we plot the total energy-loss rate, 

d\nE _ AEn(z)cr(E)c 
~dln(l + z) ~ ~E H(z) ' ^' 

as a function of energy for redshifts, 1 + z = 10, 100, and 
316. We assume AE/E ~ 1 except for the Compton- 
scattering process. Here, H(z) is the expansion rate of 
the Universe at redshift z, n(z) is the density of the target 
particle — i.e., neutral hydrogen or helium for (i) or (iii), 
free electrons for (ii) or (iv), and CMB photons for (v) 
or (vi) — and a is the corresponding cross section. At the 
high energy [E <C 13.6eV) where Compton scattering 
becomes important, a photon is not able to distinguish 
whether an electron is free or bound 2 , so we assume a 
free-electron fraction of 0.01 for processes (i), (iii), (iv), 
but treat all electrons as free for the Compton processes. 
Our result is not sensitive to the ionization fraction as 
long as the Universe is mostly neutral. We should also 
point out that in this discussion we have also neglected 
any other photon background, e.g. the infrared photon 
background which might be produced at low redshift. 



2 This was pointed out to us by Professor R. Sunyaev, see e.g. I'i.'ij 
for more detailed discussion. 
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FIG. 2: Transparency window for photons. The dark regions 
are those in which dlog E /d\og(l + z) > 1; i.e., those in which 
most of the photon energy gets absorbed by the IGM in a 
Hubble time. In the clear regions, the Universe is transparent 
to photons. 

As we can see from the Figure, high-energy photons 
(above 1 TeV at z ~ 10, 1 GeV at z ~ 300) can scat- 
ter with CMB photons or produce pairs, and they are 
largely absorbed locally, producing either an X-ray pho- 
ton which has a larger absorption cross section, or an 
electron-positron pair with high energy. 

For photons with energy around 1 GeV (and around 1 
keV at low redshift), the Universe is optically thin in the 
redshift range which we are interested in, photons lose 
most of their energy by redshifting, only a small fraction 
of the total energy is transfered to gas by scattering, and 
the scattering occurs over a wide range of redshifts. The 
details of the energy distribution depend on the injection 
energy and redshift of the photon. The formation of an 
electromagnetic shower and the resulti ng spect rum was 
investigated for both low redshift 0, 134 l35| and the 
early Universe [3(1 123 ■ The spectrum of the shower is of 
the form E~ a below the threshold energy, with < a < 2 

Fig. |21 shows the transparency window. In the dark 
regions, e?log _E/ci log (1 + z) > 1 (i.e., most of the pho- 
ton energy goes to the IGM), and in the white regions 
d log E/d log (1 +z) < 1 (i.e., the Universe is transparent 
to these photons) . The bump in the transparency window 
at (log £7, log(l + z))=(5.7, 2.16) is due to the Compton 
scattering. If the photon is injected in the transparency 
window, and remains in the transparency window as its 
energy redshifts (once a photon is injected it travels down 
and to the left on this plot), then it will free stream 
and appear in diffuse radiation backgrounds with ener- 
gies ^keV-10 TeV; otherwise, it will not appear as a dif- 
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fuse background, but will heat and ionize the IGM. Also 
keep in mind that the time interval dt oc dz(l + z)~ 5 / 2 . 
Thus, if a particle decays with lifetime longer than the 
age of the Universe, the relevant redshift for determining 
whether it appears in diffuse backgrounds is z = 0. In 
Section IV below, we will give constraints to the decay- 
particle parameter space under two extreme assumptions: 
(1) that all the photon energy goes into the IGM, and (2) 
that the photons free stream and appear as diffuse back- 
grounds. Fig. |2 must then be consulted to determine 
which limit applies for a particular decay-photon energy 
and lifetime. 

What we need to know for the ionization history is 
what fraction of the energy is converted to ionization 
energy, and how it is distributed over redshift. In the 
optically-thick case — i.e. for photons outside the trans- 
parency window — we can assume that the energy is in- 
stantly deposited. The energy deposition rate (in units 
of erg cm~ 3 s _1 ) is 



Q{z) = x l {z)nx{z)M x c 2 T 



(2) 



where Mx is the mass of the decaying particle, and 
the decay rate. If the lifetime of the particle is much 
longer than the age of the Universe and x is constant, 
then Q{z) oc (1 + z) 3 . 

In the optically-thin case, the efficiency is much lower. 
Local absorption is neglible, and a flux of high-energy 
photons is produced. These photons may interact 
with baryons by photoionization, Compton scattering, 
and pair creation, and with cosmic-radiation-background 
photons by photon-photon scattering. After interaction, 
the energy is transferred to the electrons, positrons, and 
ions which have stronger interactions with other parti- 
cles, or to photons with much smaller energy and greater 
optical depth. Here we make the following approxima- 
tion: Once a scattering happens at redshift z, the energy 
carried by the photon at that redshift is completely de- 
posited with some efficiency. This seems to be a good 
approximation, since the electrons and lower-energy pho- 
tons produced in the scattering event have much greater 
optical depths. The energy deposition rate is then 



Q(z) = Air J dE \xbnb(z)atry(E) + x-yn 7 (z)a J7 (E, z)] 
x F(E,z) (3) 

where \b, X 7 are efficiencies for converting the energy 
of that photon to ionization energy, and aby(E) and 
a 11 {E,z) are the cross sections for interacting with 
baryons and background photons, respectively, the lat- 
ter also depending on z since the background photon en- 
ergy changes. In the approximation of low optical depth 
(neglecting absorption), the flux is given by 



F(E,z) 



c 

An 



dz' 



J[^E,z>] 



(l + z')H(z') (l + z') 3 /(l + ^) 3 ' 




FIG. 3: Energy loss rate due to ionization (blue solid line) and 
inverse Compton scattering (red dashed line) of an energetic 
electron. We plot for the cases of z = 10 and z = 20. 



where J(E, z) is the emissivity at z. For simplicity, let us 
consider the case where the decay process directly pro- 
duces photons with a single energy E 1 — xMxe 2 . For 
example, if the decay products are two photons, x = 1/2. 
In other cases, the photon may not have a single energy, 
but still it is expected the energy of the photon is related 
to the mass of the decaying particle and has a narrow 
range. Then J(E,z) = N^M x c?n x (z)T X S{E - xMc 2 ), 
where iV 7 is the number of photon emitted in a decay, 
and the flux is given by 



F(E,z) 



c ( E y N^Txnxiz') 



4tt \xM x c 2 J H{z' 



1 + z l _ xM x c 
1 + z ~ E 



(5) 

Integrating over E, if the interactions with baryons dom- 
inate, then 



cn b (z)a cS 
Q( z ) Q Q ( Z >' 



(6) 



where Qo(z) is the expression given in Eq. J5J with Xi = 

NyXb, 



/dE f E \ 



(7) 



(4) 



The efficiency is roughly suppressed by a factor of t$ ~ 
nh H(z)' ' i- e - tne °Pti ca l depth for the Hubble length. If 
interaction with background radiation photons dominate, 
the suppression factor is T5 ~ n cm b(z)ca(z)/H(z). The 
redshift dependence is stronger due to the additional fac- 
tor of 1 + z in the cross section. 
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B. Electrons 

An electron can collide with and ionize atoms, or it can 
inverse- Compton scatter CMB photons, a process that 
produces an energetic photon. Those photons will be 
absorbed again, starting an electromagnetic shower. The 
energy loss dE of an electron per unit distance dx by 
ionization in a neutral hydrogen gas is given by |39j | 



dE 
dx 



mv 2 

v7 



In 



(mw V - T m ) 



1 



In 2 



2/ 2 

1 
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where 



T — 



27 m H m e v 



ra\ + m H + 2jm e iriH 



(8) 



(9) 



The energy loss by inverse Compton scattering is given 
by 



dE 4 2 

--7T = n^TCUCMBl ■ 

dt 3 



(10) 



Other forms of energy loss are relatively unimportant 
for reasonable values of parameters. For example, 
synchrotron-radiation loss is 



1 dE E ( B 

M= L05X1 ° Wfe 



(11) 



Since the energy-loss rate for inverse-Compton scattering 
is proportional to 7 2 , it dominates at high energy. At low 
energy, the ionization-loss rate is given approximately by 



dE 
dx 



2.54 x lCT ly n ff (3m7 + 20.2) eV/cm. (12) 



Since uh dc (1 + z) 3 and Uqmb dc (1 + z) 4 , inverse- 
Compton scattering dominates the energy loss at 



>20.8 



Q h h 2 \ ( 2.726 



0.022 J V Tr 
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(13) 



We plot the the energy-loss rate as a function of E 
for z — and z = 20 in Fig. [3J Generally speak- 
ing, if the electron has energy greater than ~ 100 eV 
but smaller than ~ MeV, the energy-loss mechanisms 
are collisional ionization and excitation. If the electron 
has energy greater than ~ MeV, it loses most of its en- 
ergy by inverse-Compton scattering, producing UV and 
X-ray photons. If the electron energy is E e < GeV 
or E e > 50 TeV, these photons are subsequently ab- 
sorbed by photoionization and excitation (or by pair pro- 
duction), and the decay energy is thus transferred lo- 
cally to the IGM. However, if the electron has an energy 
1 GeV < E e < 50 TeV, the up-scattered CMB photon 



has an energy in the transparency window 10 keV-10 
TeV. Thus, if the injected electron has an energy in this 
"electron transparency window," the decay energy will 
escape and appear in diffuse photon backgrounds, and 
will not transfer most of its energy to the IGM. 

For the case where the electron does heat the IGM, 
the partition of the energy among ionization, excitation, 
and heating was investigated by Shull and Van Steenberg 
[4tj|. They found that when the gas is mostly neutral, 
about 1/3 of the energy goes to ionization, about the 
same amount goes into excitation, and the rest heats the 
IGM. For a fully ionized medium, almost all of the energy 
goes into heating the gas. Therefore, we can approximate 
the fraction of energy going into ionization as, 



Xi 



Xe ~ (1 - x e )/3, Xh ~ (1 + 2x e )/3. 



(14) 



This approximation is crude but sufficiently accurate for 
our purposes. 

If the energy of the initial electron is still higher, 
E~l MeV 2 /E C MB ~ 10 10 (l + z)- 1 MeV, it can scatter 
with another photon or electron and produce an electron- 
positron pair. The electron and positron then lose en- 
ergy through inverse-Compton scattering or ionization, 
the positron eventually annihilates with another electron 
and produces 511 keV photons. To summarize, roughly 
one-third the decay energy goes locally into ionization, 
and the rest into heating the gas, unless the electron is 
injected in the transparency window GeV< E < 50 TeV, 
in which case most of the energy is carried away by up- 
scattered CMB photons. 



C. Other particles 

Protons are very penetrating particles and thus are not 
effective in transferring decay energy to the IGM. 

Other particles. Other decay products (e.g., muons, 
tau leptons, heavy quarks, gauge, or Higgs bosons) will 
generally produce showers of lower-energy particles. This 
is a complicated and model-dependent process. In the ab- 
sence of a given well-motivated candidate that decays to 
these particles, we neglect to carry out a detailed anal- 
ysis. Roughly speaking, we expect typically 10% of the 
decay to wind up in ionization energy at the decay red- 
shift, with a comparable amount going to heating the 
gas, and the rest begin carried away by neutrinos or as 
rest-mass energy of decay particles. 



III. IONIZATION HISTORY 

The decay of an unstable particle can affect both re- 
combination at z ~ 1000 and reionization at low 2, or it 
may peak at a middle redshift. However, since the first 
acoustic peak is not significantly damped, recombination 
must be rapid, and completed well before z ~ 100, when 
the Universe starts to become optically thin even at full 
ionization I27I. 
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In the presence of the decaying particle, the evolution 
of the ionization fraction x P satisfies 



dx e 1 
~~dz~ ~ (l + z)H(z) 



[R s (z)-I s (z)-I x (z)}, (15) 



where R s is the standard recombination rate, I s the 
ionization rate by standard sources, and Ix the ioniza- 
tion rate due to particle decay. This last term is re- 
lated to the energy-deposition rate Q introduced earlier: 
Ix = Q (z) I 'rib(z) I Eq , where Eq is the average ionization 
energy per baryon. At low redshift, the standard ion- 
ization sources are photons from stars or active galactic 
nuclei (AGN), I s — /*, and the standard recombination 
rate is 



R s = CmiaB{T)x 2 e n b (z), 



(16) 



where czb{T) is the case B recombination coefficient for 
gas at temperature T and density rib- Here Chii is the 
clumping factor. We take Chii = 1, appropriate for z > 

20 H3. 

The number density of decaying particles is propor- 
tional to (1 + z) 3 e~ Fxt , and its energy density is simply 
the number density times the rest mass. The particle- 
decay ionization rate is 



Ixi = e X i(z)H, with e X i{z) = Xi(z) 



M x n x (z) 
E n b {z) ' 



(17) 

To simplify the analysis, we neglect the effect of helium 
and assume jto& = m# and Eq = 13.6. The partition of 
ionization energy in hydrogen and helium depends on the 
nature of the decaying particle. The helium atom has a 
greater ionization energy and also a greater photonion- 
ization cross section, so it will probably take away more 
energy than hydrogen and produce fewer electrons. How- 
ever, it should not affect the order of magnitude of our 
estimate. The ratio rrib/Eo ~ 7 x 10 7 ; thus even only a 
tiny number of particle decays may supply enough energy 
to reionize the Universe. 

Since there are no stars present at recombination, 
CMB photons are the main source of ionization. In this 
case, a recombination to the ground state produces an 
ionizing photon which immediately ionizes another atom 
and thus produces no net change in the ionization frac- 
tion; only recombination to the n > 2 state produce net 
recombination. Assuming the number of photons in the 
2s state given by the thermal-equilibrium value, the net 
recombination rate is 



Rs — Is — C 



a B (T)xln b (z) - /3 T (1 - x e )e 



-E 2s /kT M 



where fir is the photoionization coefficient, 



C = 



1 + KAn b (l - x e ) 
1 + K{K + f3)n b {l-x e y 



(18) 



(19) 



and A = 8.23 s _1 is the two-photon decay rate of the 2s 
level. During this epoch, particle decay increases the ion- 
ization rate not only by direct ionization from the ground 



state, but also by contributing additional Lyman-alpha 
photons which boost the population at n = 2, increas- 
ing the rate of photoionization by the CMB from these 
excited states, 



Ix{z) = Ixi(z) + Ix a (z), 



(20) 



where Ixi is the ionization rate given above, and Ixa the 
ionization rate due to additional Lyman alpha photons, 



Ix a = (C-l)e a (z)H, e Xa {z) = Xa(z) 



M x n x {z)T 



x 



E a n b (z) H ' 
(21) 

Using nx — flxPc/Mx and n b = flbPc/m-b, where flx(z) 
is the fractional abundance of the decaying particle at z, 
Mx the mass of the decaying particle, and nib the mean 
baryon mass, we have 

Ixi(z) = Xi{z)^-fxT x , Ixa(z) = Xc,(z)^rfx(z)T x , 

(22) 

where fx = Qx{z)/Q,b(z). Written in this form, the 
ionization rate depends only on the fractional abundance 
of the particle and the ionization efficiency. The gas- 
temperature evolution is given by 



t dTb 8<jTanT CMB 



dz 3m e cH(z) 1 + /h c + x e 

2 K{x) 



where 



3k B H{z) 1 + /hc + x e 



nx{z) 

A(X) = XhTTlb 



(Tb 



2T h 



l CMB ) 



n b {z) x - 



(23) 



(24) 



We use a modified version of the code RECFAST [43 to 
calculate these rates and derive the ionization history of 
the gas for a given decaying particle. 

Once the ionization history is obtained, the CMB 
anisotropy can be calculated by modifying a standard 
Boltzmann code. We have used CAMB [43j for our cal- 
culation. Except for the power-spectrum normalization 
and reionization optical depth, we adopt the WMAP 
best-fit parameters for the flat ACDM model with a 
power-law spectrum i.e. {f2 m o/i 2 , Vlbah 2 , h, n s } — 
{0.14, 0.024, 0.72, 0.99}. 



A. Decaying particle with long lifetime 

The decaying particle may or may not be a major com- 
ponent of the dark matter. If T <C Hq, and the primary 
decaying particle is massive, then £l x (z) /Qb(z) is effec- 
tively constant, and the ionization history depends only 
on the energy output £ = Xifx^x- After recombination, 
the ionization fraction induced by the decaying particles 
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FIG. 4: The optical depth, IGM temperature, and ioniza- 
tion fraction as a function of redshift for standard recombi- 
nation with no reionization (black solid line) and a decaying- 
particle model with two-particle decay with Fx <C Ho and 
£ ee XifxT = 2.4 x 1CT 23 s~\r = 0.4 (red dotted line), and 
0.6 x 10" 23 s _1 ,t = 0.17 (blue dash line). 

may be estimated from the Sana equation as, 3 

x e = (e.H/aBUb) 1 / 2 oc (1 + z)^' 2 . (25) 

As the Universe expands, the physical density drops, the 
recombination rate decreases, and the ionization fraction 
increases, until at a certain point the Universe is fully 
ionized, or ionization by decaying particles is exceeded 
by stellar reionization. The contribution to the optical 
depth is then 



/odz 
—a T x e (z)n bQ (l + z) 2 oc ln(l 



+ *)• 



(26) 



In Figs. 0] we plot the ionization history for the case 
of a long-lived decaying particle (Fx <C Ho) and instant 
energy deposition, but with different energy output £. In 
this Figure, r is the Thomson optical depth between to- 
day and redshift z; Tf, is the gas temperature; and x e 
is the ionization fraction. At z ~ 1000 the ionization 
fraction drops rapidly due to rapid recombination and to 
the decrease of the ionization rate by the CMB. When 
particle decay starts to dominate the ionization rate at 
z ~ 600 — 800, the ionization fraction starts to increase 
again because now the ionization rate is constant while 



3 If there is no extra ionization, the Saha equation does not nec- 
essarily describe the ionization state as ionization reactions may 
have already frozen out because of the paucity of free electrons. 



FIG. 5: The CMB temperature and polarization power spec- 
trum 1(1 + l)Ci/(2n) for decaying particles with lifetimes 
greater than the age of the Universe. The data points with 
error bars are the binned data given by the WM AP team |44| . 
The curves are: no particle decay (black solid line), long-lived 
particle decay with £ = 2.4 x 10~ 23 s _1 (red dotted line) and 
0.6 x 10~ 23 s _1 (blue dash line). 
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FIG. 6: Same as the previous Figure, but for I < 100: 
^ = 2.4 x 10" 23 s _1 (red dotted line) and 0.6 x 10~ 23 s" 1 (blue 
dash line). We also plotted three curves for the no-particle- 
decay case (black solid line) which are almost indistinguish- 
able except for the TE polarization; from top to bottom they 
are r = 0.17, step-function reionization at z < 7, and no 
reionization. 
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the recombination rate drops. At z ~ 100, the ionization 
fraction can reach a few percent, two orders of magni- 
tude higher than in standard models. The optical depth 
increases slowly with z. 

The temperature of the IGM also starts to increase at 
z ~ 100, and continues to climb as z decreases. The rea- 
son for this is that as the neutral fraction decreases, with 
the efficiency assumed in Eq. (|14H , more and more energy 
is converted to heat at low redshift, and also at lower red- 
shift the baryons and CMB photons are kinetically de- 
coupled. The difference in the CMB and gas temperature 
produces distortions to the CMB blackbody spectrum, 
which are quantified by the Compton-y parameter, 



y = o T c 



k B (T e - T C mb) x e n b (z)dz 
m e c 2 (l + z)H(z)' 



(27) 



However, for this and all other models studied in this 
paper, we found this effect induces y < 10 -8 , well below 
the current limit [45| . 

How does the additional ionization by particle decay 
affect the CMB anisotropy? As is well known (see and 
references therein), for the temperature anisotropy there 
is only a weak effect on large scale, but on small scales 
the temperature spectrum is damped by a factor of e~ 2r . 
The division of large and small scale is determined by the 
angular size of the horizon at the reionization redshift. In 
principle, the power-spectrum normalization can also be 
determined by other measurements |46j . However, other 
parameters also affect small-scale anisotropy. To avoid 
re-fitting all the cosmological parameters, we simply fix 
all other parameters, and adjust the overall normaliza- 
tion to fit the WMAP TT and TE data. 4 

So far we have considered only energy input from par- 
ticle decay. At low redshifts, stars and quasars contribute 
a large part (if not all) of the ionizing photons. Since it is 
not the aim of this paper to provide a detailed model of 
the star-formation history, we simply illustrate the effects 
of particle decay on the CMB power spectrum by using an 
ionizing flux due only to particle decay for z > = 7; we 
then assume the Universe became suddenly and perma- 
nently reionized by standard sources below that redshift. 
We then calculate the CMB temperature and polariza- 
tion anisotropy for this ionization history. 

The renormalized CMB temperature and polarization 
power spectra are plotted in Fig. (all is) and Fig. (low 
Is). There is practically no difference in the TT spec- 
trum for the different models at high I, although there 
are small differences in the TE and EE spectrum. In 
Fig. we plot the low-/ results only; here the difference 
is more apparent. Since the overall normalization is in- 
creased, the low-/ multipoles of the high-r models are 



4 We multiply the unnormalizcd CMB power spectrum by a con- 
stant, which is then adjusted to minize \ 2 with respect to the 
first year WMAP TT (up to I = 900) and TE data (up to I = 512) 
|44| . We assume the errors are uncorrelated. 
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FIG. 7: The optical depth, IGM temperature, and ionization 
fraction for the standard no-reionization model (black solid 
line) and particle-decay-only models, all with x ~ 0-3, and 
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raised. Current observations favor low power at large an- 
gular scales m, so these models are not favored. In the 
TE data, the spectrum peaks at I ~ 10, which is again in 
contrast to the data, which is low at I ~ 10. The greatest 
difference, however, is in the EE power spectrum, which 
should easily distinguish different models. 



B. Decaying particle with short lifetime 

If the primary particle has a lifetime less than the 
age of the Universe, its density will change dramatically 
when the age of the Universe is comparable to the life- 
time. Moreover, its density today will be small. To be 
consistent with current observations of the CMB, which 
are well fit by a matter density comparable to that ob- 
tained from dynamical constraints in the present-day 
Universe, the cosmological density of the decaying parti- 
cle must be small at the time of decoupling as well. Thus, 
Qx{z) ^ Qm(z) ~ 1. In this case, 

~ %o e r x (t -t( z )) = ^xcq c -r x (t( 2 )-tc a ) ; ( 28 ) 
rib(z) flbo n bcq 

where nxeq, ^beq are the number density of decaying par- 
ticle and baryon at radiation matter equality, and the 
elapsed time is 



t(z) 



dz 



1 



1 



(1 + Z)H(Z) 3^/2(1 + ^)3/2 



(29) 
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FIG. 8: The CMB temperature and polarization power spec- 
trum. Same models as the previous Figure. 
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FIG. 9: Low-Z CMB temperature and polarization power 
spectrum. Same models as Fig. |§] The three black solid lines 
(almost indistinguishable except for the TE polarization) are, 
from top to bottom, for r = 0.17, step function reionization 
at z < 7, and no reionization. 



Unlike the long-decay-lifetime case, Fx and fx must 
now be treated as independent parameters. We show a 
few examples of the ionization history for short-lifetime 
decaying particles in Figs. and [5] The ionization rate 
is approximately constant until t{z) ~ T" 1 , after which 
it decreases rapidly. In this scenario, the ionization frac- 
tion increases slowly after recombination, just as in the 



long-lifetime case. However, because the number density 
of the particle decreases, the ionization fraction peaks 
broadly at a certain redshift and then starts to decrease 
again. The peak position depends on the lifetime. The 
models plotted in Fig. have T^ 1 = 10 14 s, 2 x 10 15 s, 
and 10 15 s, which correspond to the age of the Universe at 
z = 300, 190, and 65 respectively. Again, the ionization 
fraction can reach a few percent at z ~ 100 without jeop- 
ardizing the structure of the CMB acoustic peaks. The 
temperature of the IGM departs from the CMB temper- 
ature at redshifts of a few hundred in these cases, but 
does not increase to a very high value because of the de- 
creasing energy available for heating. The optical depth 
raises more sharply in this model, because there is more 
variation in the free-electron density at high redshift. As 
a result, the effect on the CMB is more apparent. We can 
see from Figs. [3] and |5] that the CMB temperature as well 
as polarization peaks have different shapes, especially ap- 
parent at high redshift. However, the TE correlation at 
low I is less prominent for models with short lifetime, 
and even in the model with relatively long lifetime, the 
peak is shifted to greater I, in strong disagreement with 
the WMAP result. Inclusion of low-redshift reionization 
at z — 7 results only in slight improvement. Based on 
this, particle decay with a short lifetime does not appear 
to help solve the high TE optical depth as observed by 
WMAP. 



C. Additional redshift dependence 

What if the decaying products are photons in the 
transparency window and deposit their energy differ- 
ently? As we discussed in Section II, the effect of 
photons in the transparency window in the long-decay- 
lifetime case can be described with a suppression factor of 
t,s ~ n(z)a e ec/H(z), which for baryons has a (1 + z) 3 / 2 
dependence. 

Additional dependence on the redshift may also raise 
if the density of the decaying particle does not vary as 



(1 + z) e . For example, if the decaying particle is 
relativistic, its energy density decreases as (1 + z) 4 e~ rt , 
where t' is the proper time of the moving particle. Also, 
the decaying particle could be continuously produced. 

We now consider these effects on the CMB by mul- 
tiplying £ with a factor rg(l + z) n , with ts < 1. As 
an example, we consider models with tsioq = 10~ 3 at 
1 + z = 100, and n = 0,1.5,3, and —3. Obviously, at 
least for the n = case, the ionization induced by parti- 
cle decay would be uninterestingly small if we still use the 
same parameters as in §111. A, since it is now suppressed 
by a factor of r. To see the effect of z dependence, we 
increase £ by a factor of 1000 for the n = and 1.5 mod- 
els, which cancels the small ts value we assumed. As it 
turns out, for the n = 3 and —3 models this produces 
too large a deviation during the recombination era which 
could easily be ruled out, so for the n = 3 and —3 models 
we increase £ by a factor of 100. The ionization fraction, 
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FIG. 10: The ionization fraction, temperature, and optical 
depth for models with photons in the transparency window or 
have non-standard density evolution. The curves are standard 
no reionization model (black solid line) and particle decay 
only models, all with rsioo = 10~ 3 , £ = 0.6 x 10" 24 , n = 
(red dotted line); £ = 0.6 x 10 -24 , n = 1.5 ( green short dash 
line); and £ = 0.6 x 10~ 25 , n — 3 (blue long dash line); £ = 
0.6 x 10 -25 , n = —3 (magenta dash-dotted line). 

temperature, and optical depth are plotted in Fig. 1101 

For the cases with n > 0, the additional redshift de- 
pendence (1 + z) n makes the ionization energy redshift 
away. This means that if we adopt parameters which do 
not spoil standard recombination, the effect on low red- 
shift ionization must be very small. As shown in Fig. 1111 
this does not help explain the WMAP result. 

From the above discussion, it seems that models with 
n < may produce results which arc more consistent 
with the WMAP data. For example, the n = —3 model 
tends to have relatively large effect at lower I. If, for 
example, the decaying particle is somehow associated 
with the dark energy which has a redshift dependence 
of p ~ (1 + z)°, then the n = —3 model might be real- 
ized. 

Thus, we have investigated the ionization history and 
CMB for a variety of models. It appears long-lifetime 
models may help produce the large TE polarization at 
low I. Short life time models and models with additional 
redshift dependence typically work less well, except for 
particles or other energy sources whose density decrease 
slower than ordinary matter. 

IV. OBSERVATIONAL CONSTRAINTS 

So far we have illustrated the effect of particles decay 
during the dark ages with several examples. We now 



FIG. 11: The CMB temperature and polarization power spec- 
trum for the models shown in Fig. 1101 



use measurements of the CMB power spectra and diffuse 
backgrounds to place constraints to the decay-particle 
parameter space. 

We first consider CMB constraints. We suppose that 
all of the decay energy is deposited instantly into the 
IGM, which then reionizes the Universe and changes the 
CMB power spectrum as described above. For each set 
of model parameters Tx and £ we calculate % 2 for the 
WMAP TT and TE spectrum. We have 1403 TT and 
TE data points. The \ 2 distribution for 1403 degrees of 
freedom has a a width of about 50, so the la limit is 
defined by the line Xmin + 50. The result is shown as the 
solid curve in Fig. regions above the curve are ruled 
out. Note that these results do not apply if the parti- 
cle decays to photons or electrons in the transparency 
window, as in these cases, the decay energy will not be 
deposited in the IGM but will instead propagate freely 
and appear in the diffuse radiation backgrounds. 

In this case, the £,-Tx parameter space can be con- 
strained from measurements of diffuse backgrounds z = 
0. The flux of decay radiation is then given by Eq. JSJl. 
with z = 0. We then obtain bounds to the particle den- 
sity and decay lifetime, 



N^T x n x {z) = —H(z), (30) 



where 1 + z — xMxc 2 /E. 

We now apply this result to the X-ray and 7-ray back- 
grounds. The observed cosmic X-ray background (in 
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FIG. 12: WMAP la constraints on decaying particles. Plot- 
ted are £ = xfx^x, where fx = The red solid curve 
shows constraint on the value at matter radiation equality £ cq , 
the blue dotted curve shows constraint on the value today £o • 
Note that the WMAP constraint applies if the injected photon 
or electron energy does not fall in the transparency windows 
shown in Fig. [5] and Section II. 

units of cm -2 s _1 sr _1 ) can be modeled as 

f 8 (ki/r°' 4 ' 0.2 keV < £ < 25 keV, 
F x = < 380 (^y 1 ' 6 , 25 keV < E < 350 keV, (31) 
[ 2( i J 7 )"° 7 , 350 keV < E < 2 MeV. 

The 7-ray background was measured by the the Ener- 
getic Gamma Ray Experiment Telescope (EGRET). At 
about 100 MeV, the flux is 1.45 x 10~ 5 cm" 2 s" 1 sr -1 . 
A fit for the whole energy range of 30 MeV - 100 GeV 
was H3 

F~ = 1.73 x HT 5 ( cm -2 s -1 sr -1 . 

7 ViOOMeVy 

(32) 

More recent re-analysis found a greater galactic back- 
ground, and therefore the extra-galactic background is 
lowered 0> 03- Of course, distant quasars and stars 
also contribute to this background, so here we use the 
measured background Eq. IMl't as a conservative upper 
limit. 

At each photon energy E, given the background flux 
F(E), we can derive a bound on n cq for each decay 
photon energy xMxc 2 by applying Eq. (|5U|) to the 
flux (with the restriction that the maximum redshift 
to be less than z ~ 1000). We run through the en- 
ergy range of 1 keV — 100 GeV, at each energy us- 
ing the flux limits given in Eas. (|30|l and (|32|) to de- 
rive a bound, and look for the most restrictive con- 
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FIG. 13: Constraint of £ based on diffuse X-ray and 7-ray 
background. The red solid curve shows constraint on the value 
at matter radiation equality £ eq , the blue dotted curve shows 
constraint on the value today £o- The curves are for photon 
energy (a) 100 keV, (b) 1 MeV, (c) 10 MeV, (d) 100 MeV, 
(e) 1 GeV, (f) 10 GeV, (g) 100 GeV. Note the x-ray and 7- 
ray constraints do not apply for photon and electron injection 
energies that fall outside the transparency windows. 



straint on n e q for all energies. As expected, except for 
short-lived particles, the constraint comes mainly from 
emission at z = 0. The results for decay photon in 
the transparency window xMxc 2 = 100 keV, 1 MeV, 
10 MeV, 100 MeV, 1 GeV, 10 GeV, 100 GeV, are plot- 
ted as £ eq in Fig. El O n the same Figure we also plot 
the corresponding values of today, £o- In terms of the 
decaying particle at radiation-matter equality, particles 
with short lifetimes are less constrained than the ones 
with long lifetimes. However, we do not expect any of 
these short-lifetime particles to remain today, as shown 
in Fig.HHl 

Depending on the energy of the observed photon, the 
diffuse X-ray and 7-ray background constraints are gen- 
erally more stringent than the CMB constraint except 
for short-lived particles. One must remember, however, 
that the CMB constraint applies only outside the trans- 
parency window, whereas the 7-ray constraint applies in 
the transparency window. 



V. CONCLUSIONS 

In this paper we have investigated particle decay dur- 
ing the dark ages. Such particle decay could induce par- 
tial ionization of the Universe, and thus provide a poten- 
tial alternative to early star formation as an explanation 
for the WMAP TE measurement. 
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FIG. 14: The (squared) correlation matrix for cosmologi- 
cal parameters. The parameters are the physical density of 
baryons Q.bh 2 , cold dark matter Q. c h 2 , Hubble constant h, 
power index for primordial density perturbation n a , neutrino 
density fl v h 2 , Thomson optical depth (see text) r, and the 
extra ionization parameter £24 = 10 24 f;. We assume a flat 
Universe. 



We considered how the decay energy is converted to 
ionization energy. We conclude that in many cases, a 
shower of electrons and X-ray photons are produced, in 
which case a sizable fraction (0.1 — 0.3) of the energy can 
be converted to ionization energy in situ, with compara- 
ble amount of energy going into heating the gas. How- 
ever, there are important exceptions. Photons in the 
energy range 100 keV — 1 TeV can escape, carrying with 
them most of the energy. Electrons in the energy range 
1 GeV — 50 TeV lose most of their energy by inverse- 
Compton scattering CMB photons into the above energy 
range. In these cases the ionization energy is deposited 
over a range of redshifts, and the energy deposition is 
proportional to t$ = n(z)a e sc/ H(z), where n(z) is the 
density of target particles, which can be baryons or the 
CMB photons. Because of the small value of the opti- 
cal depth ts, the decay rate must be very large to affect 
ionization. We can study reionization in these models by 
considering additional 1 + z dependence. In most cases, 
though, the models will be ruled out, as seen in Fig. 1131 
by diffuse backgrounds. 

The extra energy input from particle decays can be 
parameterized by the ionization energy input parameter 



£ = X.fx^x (£ bas the unit of s _1 ), where \ m the ef- 
ficiency fx — flx/&b, an d Tx the decay rate. If the 
lifetime of the decaying particle is longer than the age 
of the Universe, the situation is particularly simple since 
the result depends entirely on £. For short-lived parti- 
cles, one must specify both Tx and £. We studied the 
ionization history and CMB temperature and polariza- 
tion anisotropy for different cases. Although particle 
decays could partially ionize the Universe at high red- 
shift and produce a high optical depth, we found that in 
most cases they do not reproduce the WMAP result very 
well. The TE polarization does not peak at I ~ 2 but at 
I ~ 10, for example. We should pointed out however, 
this is not a unique problem with particle-decay induced 
reionization, but is also seen in other models with ex- 
tended partial reionization history. We also found that 
the EE spectrum is a sensitive probe to the ionization 
history. Furthermore, if reionization occurs at high red- 
shift, there is a change in the shape and position of the 
acoustic peaks. The ionization history is affected if the 
extra energy input has additional dependence on the red- 
shift. Typically, for an additional redshift dependence of 
(l + z) n with n > 0, the fit to CMB data is not improved, 
because the extra energy input at early times will spoil 
recombination. Models with n < may be helpful, but 
some exotic mechanism is needed for generating such a 
redshift dependence. 

We have obtained constraints on particle decays during 
the dark ages using the WMAP data as shown in Fig. ^| 
We found £ < 10 -24 s _1 for the long- lifetime case, and a 
slightly weaker bound for the short-lifetime case. How- 
ever, the short-lived particles decay at high redshift and 
we do not expect to see any left today. We also obtained 
constraints on the decaying particle from the observed 
diffuse X-ray and 7-ray backgrounds. This constraint is 
generally more stringent than the CMB constraint, but 
it actually applies to a different situation; i.e. the decay 
products are mainly photons in the energy range of the 
transparency window, where they can propagate freely 
across the Universe and contribute very little energy to 
ionization. 

The extra energy input also heats up the IGM during 
the dark age, and the temperature can rise to 10 3-4 K. 
Inverse-Compton scattering of free electrons can induce 
distortion in the CMB blackbody spectrum, but the effect 
is unobservably small (y < 10 ~ 8 ). 

If the dark-matter particle can decay, it may affect the 
estimation of cosmological parameters. To see how each 
parameter is affected, we can calculate the correlation 
matrix, which is related to the covariance matrix [50] by 



Yij — Cij/ \ZCiiCjj, 



(33) 



where the covariance matrix is given by the inverse of the 
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FIG. 15: The error ellipses for Q,t,h 2 vs. £24 and n s and £24- 
£ is positive from its physical interpretation. 



Choosing different fiducial models may affect the error 
estimates slightly. In addition to the standard parame- 
ters (physical density of baryons flbh 2 , cold dark mat- 
ter fi c /i 2 , Hubble constant h, power index for primodrial 
density perturbation n s , neutrino density f2„/i 2 , Thom- 
son optical depth r), we have added the extra ioniza- 
tion input energy £24 = 10 24 £ in the long-lived decaying- 
particle case. We will not consider the short-live case 
since it is much more model dependent. As expected, £ 
correlates strongly with the Thomson optical depth due 
to low-redshift stellar light; it will thus be difficult to 
distinguish them from CMB observations alone. Both t 
and £ correlates strongly with baryon density flbh 2 and 
the primordial spectral index n s . If a decaying particle 
exists but is neglected in the fit, then results for the val- 
ues of other cosmological parameters may be biased, and 
the error bars may be underestimated. In Fig. ^j] we plot 
error ellipses for flbh 2 — £ and n s — £ after marginalizing 
over the other parameters. 

There are other ways that particle decays during the 
cosmic dark ages could play a role in cosmology. De- 
cays might affect the recombination process j^. Parti- 
cle decays could produce a surfeit of free electrons after 
recombination; these extra electrons could then facilitate 
the formation of H2 molecules and thus potentially en- 
hance the star-formation rate. On the other hand, par- 
ticle decay may also heat up the gas, thus increase the 
Jeans mass of the primordial gas and suppress early star 
formation. The final outcome requires detailed investi- 
gation which is beyond the scope of the present paper. 
Finally if the contribution of the particle to reionization 
is significant, it may not require formation of structure 
at high redshift, thus eliminating one objection to warm 
dark-matter models. 

With future experiments, we expect to obtain more 
precise information on the ionization history than we 
have now, and should a decaying particle with t > 10 13 
sec exist, we may discover it through indirect observa- 
tions such as those discussed here. 



Fisher matrix C = F , with 

_ v [ 1 dcf T dcf T 1 dcf E dc? E ~ 

13 ~ Y [<j 2 c tt d6i 06 j + a 2 c , TE d6i 86 j \ ' 

(34) 

and Oi are the cosmological parameters to be estimated. 
We plot r 2 in Fig. HH 

In making Fig. O we have taken a fiducial model with 
the WMAP best fits with the exception r = 0.037 which 
corresponds to a sudden reionization with z re i = 6.0. 
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